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I. INTRODUCTION 

CN ■ 

I The reheating period in the Early Universe stands out as a challenge to theorists due to the close interrelationship of 
nonlinear and gravitational phenomena in its unfolding (see henceforth "papers 1-5"). The observation that, due 

to the high occupation numbers produced during preheating, most of the physics of reheating may be understood in 
^ , terms of the behavior of nonlinear classical waves has been the key to sustantial progress. The authors of papers 
—T ' 1-5 have undertaken systematic numerical simulations of the behavior of matter fields during reheating, finding that 
the reheating period is actually composed of three consecutive phases: an early one or preheating, where the dominant 
effect is the parametric amplification of matter fields out of the dynamical inflaton and gravitational backgrounds 
, an intermediate stage where the dominant phenomenon is the redistribution of energy among matter field modes 
through rescattering (in the terminology of papers 1-5), and a final stage where thermal equilibrium sets in. During 
^ I the intermediate stage, spectra of occupation numbers for the matter fields reduce to simple power laws both in the 
tJ" ■ infrared and ultraviolet limits. As noted in papers 1-5, this behavior suggests a connection between the physics of 
' reheating and the phenomena of weak turbulence |^ , but to the best of our knowledge no theoretical prediction for 
, the exponents involved is available. Our goal is to provide these theoretical estimates. 

In this paper we shall follow this same trend of ideas, by observing that, from the macroscopic point of view, 
a stochastic ensemble of classical waves may be described by a conserved energy momentum tensor subject to the 
Second Law of thermodynamics. There is therefore an equivalent fluid description, consisting of a fluid whose en- 
ergy momentum tensor and equation of state reproduce the observed ones for the microscopic fluctuations. Solving 
the dynamics of this equivalent fluid yields answers to all relevant questions concerning the behavior of the actual 
fluctuations. 

Oh, An immediate consequence of energy momentum conservation and the Second Law is that when velocities are low, 
^ ' the phenomenological fluid may be described within the Eckart theory of dissipative fluids |^ (for an analysis of the 
limitations of Eckart 's theory see [Q). It follows that it obeys a continuity equation and a curved space time Navier - 
Stokes one. The "turbulent "spectra found in numerical simulations correspond to the self-similar solutions discussed 
long ago by Chandrasekhar |11| . They appear in the discussion of decaying turbulence (which is the case relevant 
^ to cosmology), as opposed to turbulence driven by some external means. The Chandrasekhar solutions are built on 
d ' the Heisenberg closure hypothesis (see |]l3|-|l5| for a general discussion of turbulence). They were generalized to 
Friedmann - Robertson - Walker (FRW) backgrounds by Tomita et al ||l^ . These solutions agree with the Kolmogorov 
1941 theory in the inertial range [l5|] , failing to reproduce observations for very small eddies. Fortunately, we are 
most interested in the opposite limit of very large eddies, where it is trustworthy (we wish to point out that the 
applicability of Kolmogorov's spectrum to large scale turbulence should not be taken for granted ||l^). With minor 
adjustments, Tomita's analysis of turbulence decay in FRW space times also provides a solution to the evolution of 
our equivalent fluid. 

The rest of the paper is organized as follows. In next section we provide a brief summary of hydrodynamics in 
flat and expanding universes, in order to set up the language for the rest of the paper, and introduce the self similar 
solutions. In section III we proceed to discuss the equivalent fluid description of field fluctuations, and how to extract 
the particle spectrum therefrom. In section IV we place the self-similar solutions in the context of reheating. We 
state our main conclusions in the final section. We provide a rough estimate of the shear viscosity during reheating 
in the Appendix. 



II. HYDRODYNAMIC FLOWS 
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A. Flows in flat space time 



The equations governing the dynamics of a fluid in local thermodynamic equilibrium are the continuity and Navier- 
Stokes ones, which, in the case of flat space time, read: 

|^ + (U.V)p = (1) 



?^ + (U-V)U = -ivp + i/V^U (2) 
ot p 

where we have assumed incompressibility, valid when typical velocities are much smaller than the sound velocity; 
u = rj/ p is the kinematic shear viscosity. The transition from laminar to turbulent motion can be universally described 
by the dimensionless "Reynolds" number: 

(3) 

where (7 is a typical velocity and L a typical length scale. This number represents the order of magnitude of the 
ratio of the inertial to the viscous term. Low Reynolds numbers correspond to laminar motion, while high ones suggest 
turbulent behavior. 

In general, the velocity profile displays variations in space and time. This implies that the flow must be described 
probabilistically. Thus, each quantity involved in is divided in its mean value and a fluctuation from it; for 

example, we write XJ = U + u, where u stands for the fluctuating part of the velocity. In the case where motion is 
isotropic, the mean value U for the velocity must be zero, since otherwise there would be a preferred direction. 

To analyze the system's behavior, we define the two-point one-time correlation function for the velocity: 

i?y(x,x',t) = (Ui(x,t)Mj(x',t)) (4) 

In the case of homogeneous and isotropic motion, this correlation must be only a function of the time t and the 
distance between x and x', i.e. Rij{x,K',t) = Rij{r,t), where r = |x — x'| . Observe that Rii{0,t) (summation over 
repeated indices must be understood) is twice the average energy density of the fiow at time t. From (^) we obtain 
the equation that this correlation must obey, namely: 



d_ 

dt 

where 



% (r, t) = (r, t) + (r, t) + 2vV^R,j (r, t) (5) 



P^oir.t) = {p{^,t)u,{^ ,t)) - A (p(x',t)M,(x,t)) ) (6) 

p \ori OTj ' 



and 



Tij{r,t) = {ui{:x.,t)uk{:x.,t)uj{x.' ,t) - Ui{:x.,t)uk{:x.' ,t)uj{x.' ,t)) (7) 

The tensor T^- comes form the inertia term in Navier-Stokes equation and, as it involves a product of third order in 
the velocity, reflects the fact that there is not a close set of equation for the correlations of successive orders but there 
is a hierarchy of equations instead. The problem of closing that hierarchy is known as the "moment closure problem" 
[Q. Let us call $y (fc,i) the Fourier transform of Rij{r,t). Then the energy density becomes 

^RuiO,t) = J E{k,t) dk, 

where 

E{k,t)^^l <^,,{k,t) dn{-k) (8) 

is the energy density stored in eddies of size k^^ . Defining Fy as the Fourier transform of T^j, we obtain from (|^) 
the equation of balance of the energy spectrum: 
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- ^E{k, t) = T(fc, t) + 2vk^E{k, t) 



(9) 



where 

T{k,t) = ~^lTu{-k,t) dn{k) (10) 

The inertia term T{k,t) is the one that contains the mode-mode interaction, and its effect is to drain energy form 
the more energetic modes -typicaUy the bigger ones- to the ones where there is major viscous dissipation -the smaUer 
ones-. 



B. Flows in expanding universes 



For a curved space-time, in particular a Friedmann - Robertson - Walker (FRW) Universe with zero spatial curvature 
(ds^ — —dt^ + a^{t) {dx^ + dy^ -l-dz^)), the generalization of the above arguments has been considered by many 
authors |p9|-p3[. We follow Tomita et aZ.'s analysis in which they obtain the solution for the energy spectrum in 
the case of homogeneous, isotropic and incompressible turbulence. 

In a generic space time, we describe fluid flow from the energy density p, pressure p and four velocity U. The 
symmetries of the FRW solution suggest using instead the commoving three velocity = W/U'^; if W <C the 
flow is non relativistic, and if Vu =0, it is incompressible (u = (u^ , , u^)) . Later on, we shall also use the physical 
three velocity v = a{t)u. 

The corresponding continuity and Navier-Stokes equations for a Robertson- Walker background are obtained by the 
condition of conservation of the energy- momentum tensor ^ . For a non relativistic incompressible fluid, with shear 
viscosity rj = v (p + p) (but no bulk viscosity) , these reduce to: 



dp 
dt 



3-(p + p) 
a 







(11) 



dvL 

'dt 



(u-V) 



9 In ((p + p)a^ 

m 



u 



Vp 



[p + p) a? 



(12) 



where we have assumed that p + p depends only on time. For the physical three velocity v, the corresponding 
Navier-Stokes equation reads: 



9v 
'dt 



1 d\^^p + p)a^ 

a at 



Vp 



1 



a{p + p) a' 



(13) 



In obtaining (|ll|)-(|13|) we have neglected possible perturbations to the FRW metric. The corresponding equations 
considering fluctuations in the metric (g^j/ — g^^^ + /i^iy) have been obtained by Weinberg The continuity equation 
is not corrected by gravitational perturbations, while in the Navier-Stokes equation the metric fluctuations appear 
explicitly only within the shear viscosity term. It can be demonstrated that these terms involving metric fluctuations 
are negligible for scales that are inside the horizon . For scales bigger than the Hubble radius, since dissipation 
through viscosity is not effective anyway, we may still use the unperturbed Navier - Stokes equation. 

The operation of Fourier transforming in the case of a Robertson- Walker cosmology is done in terms of commoving 
wave-numbers. In doing so, the following equation for the energy spectrum is obtained: 



|i^(M) 




r(fc,t)-f 2 — 



91n ((p -I- p)a^ 
dt 



E{k,t) 



(14) 



where the relationship between E{k,t) and ^ij{k,t) as well as between T{k,t) and rij{k,t) is the same as that for a 
flat space time, if we define i?y and Ty from correlations of physical quantities, as follows: 

Rij{r,t) = (? (ui(x,i)uj(x 4- r,t)) , 



d 

Tij{r,t) = ((ui(x,t)Mfc(x,t)'Uj(x + r,t)) - {ui{x,t)uk{x + r,t)uj{x + r,t))) 



(15) 
(16) 
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C. Self similar flows in flat and expanding universes 



As we have seen in the previous section, the key element in the description of the flow is the energy spectrum 
E{k) (Eq.(^)), which is the solution of the balance equation (Eq.(||)). In it, the right hand side contains the viscous 
dissipation as well as the inertial force T{k,t. The overall effect of this term is to transfer energy from a given scale 
to smaller ones through mode - mode coupling; thus it is natural to model the action of the inertia term as a source 
of viscous dissipation, where the effective turbulent viscosity for a given mode depends on the motion of all smaller 
eddies Q. By providing closure, that is, writing this effective viscosity in terms of the spectrum itself, a closed 
evolution equation for E{k) is obtained. Concretely, Heisenberg proposed the ansatz 



T{k',t) dk' = 2v{k,t) I E{k',t) k' 2 dk' (17) 



where 



^(^k^t) = Aflat £ \j^^dk' (18) 
and Aflat is a dimensionless constant. With this hypothesis (known as the Heisenberg hypothesis) as the solution to 



the closure problem, Chandrasekhar 
there is a stage in the decay where t 



11| has obtained the energy spectrum for decaying turbulence, assuming that 
le bigger eddies have sufficient amount of energy to maintain an equilibrium 
distribution, thus requiring that the solution for the spectrum should be self similar. With this consideration into 
account he obtained an energy spectrum: 



E{k.t)^^r^,MF(p^] (19) 



A)iat^tl\! t \koVk 

where fco and Iq are initial conditions (namely, the wave number corresponding to the bigger eddy and its typical 
time of evolution). The function F obeys the equation 



, „ ,dF(x') 
F{x') -x' ^ ' 



dx' 




dx' = I vklto + / ^ ,Y;' dx' / F{x')x' ^dx' (20) 



which predicts a Kolmogorov type behavior for an inviscid fluid {R oo, R — (yJiQio) ^) in the ultraviolet limit: 

F{x) const x-^/^ (i^ = , x-^oo) (21) 

While for nonzero viscosity: 

F{x) const x-"^ {v^Q , x ^ oo) (22) 

In the infrared limit, F has the universal behavior F{x) = 4x (x <C 1), and thus we find a linear energy spectrum 
for k^/t < fcoV^o- 

Chandrasekhar's self similar solutions are easily generalized to flows in expanding Universes. The dependence on 
time and wave-number for the self similar energy spectrum is |l6| 

E{k,t) = ]^vUt)mF{m (23) 
where A and vt are respectively the Taylor's microscale and an average turbulent velocity, defined as: 

The second equation implies the normalization condition 

F(x')dx' = 1 (25) 



4 



To obtain a self similar flow reducing to Eq. (n9^ in the flat space limit, we must require the time evolution laws: 



A2(i) = + 10 



dt vt = vti 



The equation which determines de function F{Xk) in ( |23| ) turns out to be 

,dF{x') 



F{x') - x' 



dx' 



dx 




-J^dx' I f F{x')x' ^dx' 



(26) 



(27) 



where A is a constant. This equation has the same structure as in flat space time, Eq. (pOD, which means that 
assuming Heisenberg's hypothesis the spectrum is linear in k for length scales much bigger than the Taylor's microscale. 

Observe that for flat space time, the proportionality between the integral up to a certain wave number k of the 
inertia and the viscous forces is given by (^7|) and (|l8|). In the case of a FRW space time, the autosimilar solution 
required by Tomita et al. ( p3| ) needs a time dependent Vcurv in eq. (|l7|), defined by the analog of (18) with A fiat 
replaced by Acurv — 5Ar/a^. 



D. Solving for the spectrum 



Let us analyze in more detail the solutions of eq. (p7|). We assume the normalization eq. 
x — *■ oo limit in eq. ( p7|) we find 



where 



/>oo 

/ F{x')x'^dx' = 5 

^0 



Taking a derivative of eq. (g^) we get 



l-x^^x'U + AiG^H] 



S]). By taking the 
(28) 

(29) 

(30) 



H = 



F{x')x' ^dx' 



(31) 



Let us consider first the a; ^ limit. Assume F cx x". The left hand side of eq. ( ^91) tends to a finite limit 1 — a. 
On the right hand side, G and H behave as x'""^^/^, so if a > 0, this side goes to zero. We must therefore have 
a = 1, and 



F^^Cx, x -> 



(32) 



were C is some constant. 

In the a; — > oo limit, assume again a power law behavior F cx x~^ . Now G — > 0, so we must have H —> 2/5A. From 
eq. (p8|) we know that in this limit H ^ 5/ V Fx"^ , so we must have /3 — 7 and 



F . 



/25^y 



y 2 J ' 

Taking into account both limiting behaviors and eq. (^5|), the function F may be approximated as 



F\x\ = 



[a + [3x^] 



25^7r2 
128 



1/3 



(33) 



(34) 



5 



III. EQUIVALENT FLUID FOR FIELD FLUCTUATIONS 



After establishing the basic necessary notions for the description of hydrodynamic flows, our goal is to associate an 
equivalent fluid description to field fluctuations, and to derive the particle spectrum therefrom. Our flrst step is to 
obtain the energy density, pressure and velocity of this fluid as functionals of the quantum state of the field. 

For simplicity, we shall consider the theory of a single, self-interacting scalar field 4>, minimally coupled to gravity. 
The action is 

S = [ d'xV^K^] - ^ M I (35) 



_ 2 ^ 

where V [(/)] ia renormalized effective potential. The energy-momentum tensor is associated to the Heisenberg 
operator 

T^'' - = a^<^a''</. - g^'- 1 dpct^d"^ + y (36) 

The macroscopic behavior of the field, however, may be described in terms of a c-number energy-momentum tensor 

Tf"" = Tg" + T^" (37) 

where 



'c 



and Tg is a stochastic component with zero mean and self-correlation 

iTrn'')s = \ ({^r'^r}>Q - ^r^^r (39) 

In these equations, ()g denotes a stochastic average, while ()q is the average with respect to the quantum state 
of the field. Following Landau, we define the fluid four-velocity and energy density p as the (only) time-like 
eigenvector of T^^ and (minus) its corresponding eigenvalue 

T^'-'C/i, = -pU^ (40) 

Introducing the pressure p = p (p) as given by the equilibrium equation of state (our theory does not have a 
conserved particle number current, and therefore the equation of state is barotropic), we may decompose 

Tf"" = pU^'U•' + pA^"^ + T^"' (41) 

where A'"^ = (7'^" ~\-U^U'^ and by construction t'^'^Ui, = 0. Since r'^'^ vanishes by definition in the equilibrium state, 
it may be parametrized in terms of deviations from equilibrium. Remaining within the so-called first order formalism 
pa,E6|, we may write 



r'"' = -riH'''' - CUPpAf""] ??, C > (42) 

where 



2 



2 



(43) 



and 77 and C are the shear and bulk viscosity coefficients, respectively. 

Let us decompose each quantity in a mean component (denoted by a C subscript) and a fluctuation (denoted by 
a. S). If the quantum state shares the symmetries of the FRW background, then = 0. Since — —1 holds 
identically (as oppossed to " in the mean" ) we must have 



U'cT + mT)-a^Um)^^l (44) 



2U°cU°s - a' [WsUh - {UhUh)s] = (45) 
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The second equation shows that Ug is a higher order fluctuation with respect to Ug. If we remain at hnear order, 
then, we may approximate = = 1. Observing that all mean values are homogeneous and isotropic, we see that 
r"* is also a higher order fluctuation. We find 

T« = r« = (p + p)^C/^ (46) 

and therefore the velocity correlation 

"■'<''''^i(;fi:w^{''«'''''''2''°-''})« '''' 

This is the key equation linking the quantum and stochastic descriptions. To estimate the velocity correlation, let 
us assume that, after integrating out the hard modes, the soft modes of interest may be described in terms of quasi 
free, long lived excitations with an effective mass (t) . Then 



-1 



The fluctuations are Gaussian to a very good approximation, and therefore 



dt^d,(f> (48) 



{r« (f,t) ,r°^ (0,<)})^ = (^-jL.) {d^^,G+dl,G+ + dl,G+dl,G+ + (f^ -r)} 



(49) 



where (a;,x') is the positive frequency propagator 

G+ ((r, t) , (0, t)) = (0 (r, t) (0, t))^ , (50) 

{di^t stand for derivatives with respect to the first argument, while dii ^t' stand for derivatives with respect to the 
second argument of G+). Let us decompose the soft field into modes 

^e^'%W (51) 



At any time to we may introduce positive frequency adiabatic modes defined by 

1 / . r 

— , cxp < —t 

V^^k (t) I Jto 

where 



fk (t) = _L^ cxp ( -^ / dt' LUk it') )■ (52) 



coUt)^^^+M^t) (53) 

and decompose the mode amplitude 0j (t) into positive and negative frequency components 

(bj: (t) = ,h (t) Aj: (t) + A! (t) (t) (54) 

a*</.£ (i) = -tiuk (t) [fk it) Aj: (t) - /,! (t) Al^ it) } (55) 

Let us define the spectrum 

n,it) = {Alit)A~,it))^ (56) 
(because the quantum state is isotropic, the spectrum depends only on k) and assume that 



- (57) 
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(this happens, for example, if the different modes acquire random phases through interaction with an environment). 
Then 

j-^ e''^ {fk it) ft it') (1 + n, {t)) + ft {t) fu {t') nu {t)} (58) 

And 

—3 e^'^^^ {-zu, [t)) {fu {t) ft [t') (1 + it)) - ft (t) A. it') n, {i)} (59) 

(zTTj 

Observe that only the vacuum part contributes in the coincidence limit t' t. In the large occupation numbers 
regime we are interested in, this is negligible, and we get 

({T«(f,t),r°^"(0,i)})^^ (^^) {df^,G+d?,,G+ + {r^-r)} (60) 



where 



€'G+ - I N/c {t) (61) 



r^'j {2^f Wit) 



(Pk 
(2^ 



e^^^^Uk {t)nk it) 



We may now write down the Fourier transform of the velocity self correlation 



(fc,<) = 



3a^t)ip + p)'c{t)J i2n)' W 



(Pp UJp (t) 



k — p 



(62) 



(63) 



In principle, this is an integral equation relating the spectrum to the energy self-correlation. It may be simplified 
as follows. For small p < k, the integral reads 



V^fc (t) 



Uk (t) 



(Pp 
(27r 



-3 UJp (t) Up (t) 



while for large p we get 



(t) 



(64) 



(65) 



Since the spectra we are considering fall off much faster than p ^ (or even p ^) at large p, we may estimate that the 
contribution from cq. (|6^) is larger than eq. (p5|); moreover, extending the integral to the whole momentum space, 



d^p 
{2nf 



UJp (t) Hp (t) p (t) 



and we get 



(fc,t) 



Pit)S'' 



3aHt)ip + p)l{t) Wit) 



Uk it) 



Finally, we may relate the particle spectrum to the turbulent energy spectrum 



E{k,t) 



2ttp [t) 



aHt){p + p)l{t) Wit) 



Hk it) 



(66) 



(67) 



(68) 
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IV. REHEATING 



In the previous sections, we have analyzed on one hand self similar turbulent flows in expanding universes, and on 
the other hand have given the rule to translate the turbulent energy spectrum into a particle number distribution. 
We must now show that the foregoing analysis is relevant to plausible models of the reheating period, and use it to 
predict the likely shape of the final particle spectrum. 

At this point, it is convenient to be more precise about the model of reheating we have in mind (for further details, 
see papers 1-5). We assume that inflation is driven by an inflaton field <i> with an effective potential which may be 
parameterized as V ($) ~ Ai„/$''. The self coupling Xmf ~ 10~^^,and at the end of inflation <i> ^ irip ~ lO^^GeV (in 
natural units). This is the dominant contribution to the energy density, so it fixes the scale for the Hubble constant 
H: mlH^ ~ \/($). 

During preheating, a large fraction (nearly all) of this energy is transferred to the matter fields. These are described 
by an effective degree of freedom which self interacts with an effective dimensionless coupling constant g. For 
simplicity, we shall model this self interaction as a gcf)"^ theory. At times to at the end of preheating, matter excitations 
are distributed with a smooth spectrum Uk ~ Nf (fc/fco), where fc is a comoving wave number and the characteristic 
momentum H (meaning that the relevant modes are inside the horizon). We wish to find the form of the 

spectrum as follows from the assumption that its further evolution will be self-similar. 

'4 ~ ^inf^ 



Let us estimate the energy density in the matter fields at as po ~ iVfcg ~ Ai„/$^. The self interaction of the 



matter fields induces a mass M^. We shall assume > Af^, in which case 

— Uk ^ gNkl ^ g {po/kl) . (69) 

We shall assume that this is the dominant contribution to the matter self-energy. The condition fcg > is 
equivalent to 

gN < 1 (70) 
(this equation will be important below). Observe that 



M'^ gPo f 9 \fPo\fV{^)\ f g \fH\ (m 



k^ /Cq y ^inf J \ ^0 / \ ) V Airi/ 



Since i?/fco is already less than one, this is not unduly restrictive on g (see papers 1-5). Also observe that 
redshifts with the cosmological evolution, so these estimates are not affected by Hubble fiow. 
Integrating equation ( |68| ) over /c, we find the mean velocity in the equivalent fluid flow as 

4 - < (72) 

(Po +Po) 

where 

'' IP- J_ /lf2 ^ ' 

{kophys — ko/a (t)) is the speed of sound, so the fiow may be regarded as incompressible. The shear viscosity of the 
matter fields may be estimated as (see appendix) 

k^ 

- - (74) 



leading to a kinematic viscosity 



Ng 



— r- (75) 

P + P {gNfko 



With a typical velocity of order one, and a typical lenght of order fcg ^, we get the Reynolds number R ^ (gN)"^ < 1 
(cfr. eq. dT^)). 
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Of course, we do not have a microscopic justification for Heisenberg's closure, so we shall take A in eq. ( |27| ) as a 
free parameter. As it turns out, agreement with the numerical results reported by papers 1-5 is obtained for A ^ 3, 
which, as expected, corresponds to Reynolds numbers of order of one. 

It's worth to point out that Chandrasekhar's solution with the function F obeying eq. (|27| ) is an exact solution 
that relies only on Heisenberg's closure hypotesis eq. (JT?!). So if we trust this hypotesis, we can get the spectrum by 
solving eq. (|2^) for any Reynolds number. 

The distribution of occupation numbers is found from eg. (|68|), where the energy spectrum is given by eq. (p3|). In 
the light field limit Af^ <C fcg, we find cuk ~ k. Then eq. (|3^ ) implies Uk ~ for k ^ 0, and eq. ( ^3| ) implies A:^^" 
for fc — > oo. This theoretical prediction for the exponents involved is the main result of this paper. 

In Fig. 1 we show the full particle spectrum based on the approximation eq. (^J) for the function F. We have 
scaled the plot to make it easiest to compare with papers 1-5. Momentum is measured in units of a^^lO^^GeV (where 
a is the scale factor) and we have chosen A ~ a 1.8 x 10~^^GeV~^ . Observing that in the ultrarelativistic limit 
the speed of flow vt ^ 1 and p ^ p, integrating the particle density times shows that the total energy density in the 
flow is p ^ 10^ X (lO^^Gey) . The equivalent black-body temperature is then somewhat less than lO^^GeV, which is 
a reasonable value for the reheating era, and high enough to justify the neglect of all masses. 

The analysis leading to the Chandrasekhar solutions begins from the Navier - Stokes equations for the fluid, which 
are equivalent in this context to energy - momentum conservation. Therefore our model requires that it be possible to 
assign to the fluid an independently conserved energy - momentum tensor. This is not exactly the same as requiring 
that the homogeneous mode has totally decayed, but it does mean that there is no significative energy transfer from 
the homogeneous mode to the fluid, either through parametric amplification or otherwise. The plots of total particle 
density and effective masses presented in ref. [Q suggest that this condition is met early, then follows an intermediate 
stage dominated by rescattering, and finally the thermalization stage. 

Comparison with the results in papers 1-5 is meaningful only in the intermediate phase. Decaying turbulence is 
necessarily a transient phenomenon. As time evolves, we expect the field will eventually thermalize, and the spectrum 
will get closer to a Rayleigh-Jeans tail, ~ k~^ when masses are negligible. The several plots presented in ref. 
Q, where indexes go from 1.7 to slightly over 1, capture the transition from turbulence to equilibrium. Since the 
same plots show that earlier spectra are steeper (see also Fig. 1 in ref. |^) this is in satisfactory agreement with the 
prediction from self-similar flows. 
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Fig. 1: Log-Log plot of the particle spectrum, as given by eq. (|68|). The energy spectrum is given by (|23|), where 
the function F is given by eq. (jsj). We have chosen A = 3, the Taylor microscale A ^ 0.18 x 10~^^GeV~^ , and have 
scaled the spectrum to make it easiest to compare against the results presented in |^ and 



V. FINAL REMARKS 



In this paper, we have shown that the self similar flows studied by Heisenberg, Chandrasekhar and Tomita may be 
used to provide an interpretation of the "turbulent "spectra found in papers 1-5. The hydrodynamic model predicts 
scale invariant spectra Uk ^ fc~" both in the infrared and ultraviolet limits, with a ~ 2 in the former, and 10 in the 
latter regime. Agreement with the early time results presented in papers 1-5 is satisfactory. 
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The connection of hydrodynamics to the behavior of fluctuations during reheating has interest of its own, as it 
provides an alternative to brute force quantum field theoretic calculations, and also yields physical insight on the 
macroscopic behavior of quantum fields in the Early universe. The equivalent fluid method may be used to advantage 
also in other regimes, such as the inflationary period itself ||29| ]. Moreover, it opens up a wealth of new phenomena, 
such as intcrmittence [ p^ and shocks |]30[|, which are not apparent in the customary treatments. We will continue 
our research in this field, which promises a most rewarding dialogue between cosmology, astrophysics, and nonlinear 
physics at large. 
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VII. APPENDIX: ESTIMATE FOR THE MATTER VISCOSITY 

In order to estimate the viscosity for the matter fields we assume that the correspondence eq. ( p8| ) between the fluid 
and particle spectra allows us to associate a solution of the hydrodynamical evolution equations for the former to a 
self - similar solution of the kinetic equation for the latter (cfr . [pi ). We may then estimate the transport coefficients 
by adopting the same methods usually applied in equilibrium mW. These are discussed in detail, in the quantum field 
theory context, in references [ p8[ . 

Let n1 ~ Nf {X), X = W^kf^/ko be a solution to the covariant Boltzmann equation 

fc^V^rifc = Icoi [uk] (76) 

Since we are interested in the high energy regime where shear viscosity is much larger than bulk viscosity, we may 
use the Boltzmann collision integral 

Icoi [nki] 9^ J Dk2Dk3Dk4 5 (fci + k2 - k^ - k^) 

{(1 + nki) (1 + nk2) nkzUki - (1 + ^^fe3) (1 + "m) nkink2} (77) 

where Dk ~ d^fc S (k^)- To estimate the shear viscosity we assume the four-velocity is slightly inhomogeneous, 
and seek a solution 

nk-nl + nl (78) 
where the new term satisfies the equation (which we render only schematically) 



By simple power counting, we estimate 



Sn 
And so 



[nl] ^ g^N^kl (80) 



1 df 



The new term induces a correction to the energy - momentum tensor 

= I Dk hkjTil (82) 
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Equating this to rjUij and repeating our power counting analysis, we obtain 



l3 



(83) 



as in eq. (74). In equilibrium, ko ^ N ^ I, and we obtain the same result as in refs. 
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